Method for indicating fatigue damage of a metal object

ABSTRACT

A method for indicating fatigue damage rate of a hardened metal object in relation to load cycles, N, exerted on the hardened metal object, wherein the hardened metal object presents a temperature essentially corresponding to the operating conditions of the hardened metal object is disclosed. It comprises calculating the fatigue rate based on an effective activation energy parameter for the dislocation climb process, Q, shear stress amplitude, T, the absolute local temperature of the hardened metal object, T, and load frequency, f. Methods for indicating fatigue and predicting life of a metal object are disclosed.

TECHNICAL FIELD

The invention concerns methods of indicating fatigue damage and damage rate of a metal object, especially a hardened metal object, and in particular a hardened metal bearing component.

BACKGROUND

Fatigue, one of the most common material degradation mechanisms in industry in general and in the bearing industry in particular, occurs when material experiences lengthy periods under repeated or cyclic stresses which can lead to failure at stress levels much lower than the tensile or yield strength. It has long been recognised that nearly 90% of industrial component failure takes place due to fatigue. Hence the importance of evaluation of fatigue damage in metallic components and building a solid understanding of the fatigue phenomenon, aiming at preventing fatigue failures from occurring.

One example in this direction is EP1184813A2. EP1184813A2 describes how a dynamic equivalent load P is calculated from data information of a rolling bearing. Thereafter, a reliability coefficient a_(l) is determined, a lubrication parameter a_(L) corresponding to a used lubricant is calculated, and a contamination degree coefficient a_(c) is determined in consideration of a material coefficient. A fatigue limit load P_(u) is calculated on the basis of the data information. Thereafter, a load parameter {(P−P_(u))/C}·1/a_(o) is calculated. On the basis of the lubrication parameter and the load parameter {(P−P_(u))/C}·1/a_(o), a life correction coefficient a_(NSK) is calculated with reference to a life correction coefficient calculation map. The bearing life L_(A) is calculated as L_(A)=a_(l)·a_(NSK)·(C/P)^(P).

Unfortunately life models commonly approach the fatigue problem in a simplified manner and commonly treat the metal in a highly idealised way. The models therefore exhibit limitations and there seems still to be room for improvements as regards the understanding of the characteristics of metals, especially as a bearing material, the life models, and the influences from the operating environment on the characteristics of bearing material.

SUMMARY

According to one aspect, the present invention relates to a method for indicating fatigue damage rate of a hardened metal object and in particular a hardened metal bearing component.

According to another aspect, the present invention relates to a method for indicating fatigue damage of a hardened metal object and in particular a hardened metal bearing component.

According to still another aspect, the present invention relates to a method for indicating at least one of an effective activation energy parameter for the dislocation climb process, shear stress amplitude, the absolute temperature of the hardened metal object, load frequency, and number of load cycles exerted on the hardened metal object.

According to a further aspect, the present invention relates to a computer program product implementing at least one of the previous aspects.

According to still a further aspect, the present invention relates to a computer readable medium containing program instructions according to any of the previous aspects.

The present invention has its basis in a metal physics description of how hardened metal materials behave under high cycle fatigue conditions, including not only the effects of the Hertzian contact stress field, but also the effects of the operating temperature, superimposed (hoop and residual) stresses and speed. The term fatigue includes at least one of rolling contact fatigue (RCF) and structural fatigue, such as rotating bending fatigue, torsion fatigue, uniaxial fatigue, including push-pull fatigue, and multiaxial fatigue. The present invention deals with plastic damage rate, accumulated fatigue damage, cumulative fatigue, strain damage, and cumulative strain damage.

According to an aspect of the invention a method for indicating fatigue damage rate of a hardened metal object in relation to load cycles, N, exerted on the hardened metal object, wherein the hardened metal object presents a temperature essentially corresponding to the operating conditions of the hardened metal object is disclosed. It comprises calculating the fatigue rate based on an effective activation energy parameter for the dislocation climb process, Q, shear stress amplitude, T, the absolute local temperature of the hardened metal object, T, and load frequency, f. Further also methods for indicating fatigue and predicting life of a metal object are disclosed

According to one aspect of the present invention, a method for indicating fatigue damage rate of a hardened metal object in relation to load cycles, N, exerted on the hardened metal object will now be disclosed. The hardened metal object has a local temperature. Local temperature is here defined as a temperature essentially corresponding to the operating conditions of the hardened metal object. The method comprises calculating the fatigue damage rate on a load cycle base dγ/dN or a time base by using an alternative formalism, which will be further elaborated below. The calculating is based on an effective activation energy parameter for the dislocation climb process, Q, shear stress amplitude, T (expressed in Pa), the absolute local temperature of the hardened metal object, T, and load frequency or rotational speed, f. In an embodiment the load frequency is constant. In another embodiment, the load frequency varies over time. In an embodiment, the fatigue damage rate is indicated to a user. Each one of the shear stress amplitude and the temperature are constant or variable over time.

In an embodiment, and in case of the shear stress and/or temperature is constant or non-constant, the fatigue damage rate may be expressed as dγ/dt=A<T−T _(u)>^(c) e^((−kQ/RT)), where

-   -   A is A₀ (H_(Vref)/H_(V))^(d), where A₀ is a constant larger than         zero but smaller than 10⁻⁸, H_(Vref) is the Vickers hardness for         a reference metal, H_(V) is the Vickers hardness for the         hardened metal object,     -   <T−T _(u)> is zero for T≦T _(u) and T−T _(u) for T>T _(u), where         T _(u) is the fatigue limit for the hardened metal object,     -   c is a constant in the interval 6 to 22,     -   d is a constant in the interval 6 to 22,     -   k is a constant between 0.50 and 1.50,     -   T is the absolute temperature of the hardened metal object at         the point of contact, and     -   R is the universal gas constant.

However, substituting t=N/f in the expression above (and thus dt=(1/f)dN), an embodiment is reached, in which the step of calculating comprises calculating dγ/dN according to dγ/dN=A< T−T _(u)>^(c) e ^((−kQ/RT)) /f, where

-   -   A is A₀(H_(Vref)/H_(V))^(d), where A₀ is a constant larger than         zero but smaller than 10⁻⁸, H_(Vref) is the Vickers hardness for         a reference metal, H_(V) is the Vickers hardness for the         hardened metal object     -   <T−T _(u)> is zero for T≦T _(u) and T−T _(u) for T>T _(u), where         T _(u) is the fatigue limit for the hardened metal object,     -   c is a constant in the interval 6 to 22,     -   d is a constant in the interval 6 to 22,     -   k is a constant between 0.50 and 1.50,     -   T is the absolute temperature of the hardened metal object at         the point of contact, and     -   R is the universal gas constant.

In an embodiment, the absolute temperature T expressed in degrees Celsius, t, is in the interval of −55 to 250 degrees Celsius, where t=T−To and To=273.15 degrees Kelvin. In other embodiments the temperature, t, is in the interval of zero to 250 degrees Celsius, in the interval of 15 to 200 degrees Celsius, in the interval of 30 to 150 degrees Celsius, in the interval of 50 to 120 degrees Celsius, in the interval 65 to 110 degrees Celsius, or in the interval 75 to 100 degrees Celsius. In some embodiments, the absolute temperature, T, of the hardened metal object is below 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.

In some embodiments, the method further comprises calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q₀−ΔVσ, where ΔV is a material constant, and σ is the normal stress field tensor. According to the invention it has also been determined that the fatigue damage process is also influenced by the superimposed (static or dynamic) normal stress field. In other embodiments, σ is the hydrostatic stress, which in an embodiment, is calculated according to: σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are the normal stress components.

In another aspect of the invention, integrating the fatigue damage rate relation according to the first aspect over time, or more specifically the variable of number of load cycles, using the actual operating conditions, gives the cumulative fatigue damage, i.e. the integral of the first aspect. A method for indicating cumulative fatigue damage, γ, of a hardened metal object at a temperature essentially corresponding to the local operating conditions of the hardened metal object and under constant stress amplitude will now be disclosed. It comprises calculating the cumulative fatigue damage, γ, based on effective activation energy parameter for the dislocation climb process, Q, shear stress amplitude, T, the absolute local temperature of the hardened metal object, T, load frequency, or rotational speed, f, and number of load cycles, N, exerted on the hardened metal object.

In a case of the shear stress and temperature being constant, the fatigue damage may be expressed as γ=A<T−T _(u)>^(c)e^((−kQ/RT))t. This is the integral of the expression concerning the time based fatigue damage rate discussed above.

However, using also, for a constant load frequency, the substitution t=N/f, an embodiment is reached, in which the step of calculating comprises calculating γ according to: γ=A> T−T _(u)>^(c) e ^((−kQ/RT)) N/f, where

-   -   A is A₀(H_(Vref)/H_(V))^(d), where A₀ is a constant larger than         zero but smaller than 10⁻⁸, H_(Vref) is the Vickers hardness for         a reference metal, H_(V) is the Vickers hardness for the         hardened metal object,     -   <T−T _(u)> is zero for T≦T _(u) and T−T _(u) for T>T _(u), where         T _(u) is the fatigue limit for the hardened metal object,     -   c is a constant in the interval 6 to 22,     -   d is a constant in the interval 6 to 22,     -   k is a constant between 0.50 and 1.50,     -   T is the absolute temperature of the hardened metal object at         the point of contact, and     -   R is the universal gas constant.

In some embodiments, the absolute temperature T expressed in degrees Celsius, t, is in the interval of −55 to 250 degrees Celsius, in the interval of zero to 250 degrees Celsius, in the interval of 15 to 200 degrees Celsius, in the interval of 30 to 150 degrees Celsius, in the interval of 50 to 120 degrees Celsius, in the interval 65 to 110 degrees Celsius, or in the interval 75 to 100 degrees Celsius, where t=T−To and To=273.15 degrees Kelvin. In other embodiments, the temperature, T, is below 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.

In some embodiments, the method further comprises calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q₀−ΔVσ, where ΔV is a material constant, and σ is the normal stress field tensor. This is due to that the fatigue damage process is also influenced by the superimposed (static or dynamic) normal stress field. In some embodiments, σ is the hydrostatic stress, and in other embodiments, σ is calculated according to σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are the normal stress components.

Now turning to still another aspect and assuming that metal failure occurs when the highest stressed point has accumulated a (non-specified) critical degree of plastic damage, then the formula presented in the second aspect may be used to predict the metal fatigue life. According to this aspect, a method for indicating at least one of an effective activation energy parameter for the dislocation climb process, Q, shear stress amplitude, T, the absolute local temperature of the hardened metal object, T, load frequency or rotational speed, f, and number of load cycles, N, exerted on the hardened metal object, at a temperature essentially corresponding to the operating conditions of the hardened metal object and under constant stress amplitude is disclosed. In an embodiment, the method may also be used to indicate normal stress field or hydrostatic stress. The method is based on the following relation N=Cf/(<T−T _(u)>^(c) e ^((−kQ/RT))), where

-   -   C=C₀ (H_(V)/H_(Vref))^(d), where C₀ is a calibration constant         representing operating conditions for a reference metal object         having a known fatigue life in terms of estimated maximum number         of load cycles, H_(Vref) is the Vickers hardness for a reference         metal and H_(V) is the Vickers hardness for the hardened metal         object,     -   <T−T _(u)> is zero for T≦T _(u) and T−T _(u) for T>T _(u), where         T _(u) is the fatigue limit for the hardened metal object,     -   c is a constant in the interval 6 to 22,     -   d is a constant in the interval 6 to 22,     -   k is a constant between 0.50 and 1.50,     -   T is the absolute temperature of the hardened metal object at         the point of contact, and     -   R is the universal gas constant.

Thus, C₀ is a constant that represents the reference operating conditions for the metal object having a known life (for example from endurance testing) and which was used to calibrate the model.

The industrial applicability of the formula above is of course not limited to estimating N, but, for instance when a prediction of N may be known before hand, to solve the equation in view of for instance the local temperature.

In some embodiments, the absolute temperature T expressed in degrees Celsius, t, is in the interval of −55 to 250 degrees Celsius, in the interval of zero to 250 degrees Celsius, in the interval of 15 to 200 degrees Celsius, in the interval of 30 to 150 degrees Celsius, in the interval of 50 to 120 degrees Celsius, in the interval 65 to 110 degrees Celsius, or in the interval 75 to 110 degrees Celsius, where t=T−To and To=273.15 degrees Kelvin. In other embodiments, the temperature, T, is below 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.

In some embodiments, the method further comprises calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q₀−ΔVσ, where ΔV is a material constant, and σ is the normal stress field.

This is due to that the fatigue damage process is also influenced by the superimposed (static or dynamic) normal stress field.

In some embodiments, σ is the hydrostatic stress, and in other embodiments, it is calculated according to σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are the normal stress components.

According to a further aspect, a computer program product is disclosed. It is loadable into the internal memory of a computer, comprising software code portions for performing step(s) of any of the previous aspects, when run on a computer.

According to still a further aspect, a computer readable medium containing program instructions for execution on a computer system, which when executed by the computer system, cause the computer system to perform step(s) of any of the previous aspects.

In a non-limiting embodiment the aspects of the present invention may be utilised in remote control, or remote monitoring, of machinery comprising hardened materials by measuring and monitoring the variables discussed above.

The different enhancements of the invention as described above can be combined in any desired manner as long as no conflicting enhancements/features/characteristics are combined.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention will now be described in more detail for explanatory, and in no sense limiting, purposes, with reference to the following figures, in which

FIG. 1 shows a schematic illustration of a device suitable for executing the methods according to present invention,

FIG. 2 shows a schematic illustration of a fatigue test rig,

FIG. 3-5 shows diagrams indicating the relation between fatigue damage, in terms of surface groove depth, number of load cycles, and the absolute temperature of the hardened metal object.

DETAILED DESCRIPTION

In order to clarify the method and device according to the invention, some examples of its use will now be described in connection with FIGS. 1 to 5.

High cycle fatigue response can be seen as being a function of the applied stress amplitude and the number of stress cycles, i.e. a mechanical approach. The fatigue life is then to some extent seen as a probabilistic process that can be treated using statistical methods. Correction factors for any applied mean stress, lubricant film thickness and lubricant contamination are then applied in metal life models for rolling contact fatigue (RCF). The relationship between the applied bearing load (P) and the life (L), expressed in million revolutions, has the form: L=a _(SLF)(C/P)^(p)≈Constant/<T−T _(u)>^(c), where a_(SLF) is the stress life factor, C is the basic dynamic load rating, p is the load ratio exponent, T is the subsurface shear stress amplitude, T _(u) is a fatigue limit shear stress and c is the shear stress exponent. Here <T−T _(u)> is zero for T≦T _(u) and T−T _(u) for T>T _(u). However, it should be noted that this does not treat fatigue damage as being a microstructural effect and therefore only describes the material rather as being a “continuum” having predefined (and non-changing) properties. Still metallographic techniques, including microscopy and X-ray diffraction methods can be used to monitor the accumulation of fatigue damage in mechanical objects, such as bearings. This is an implicit use of a changing microstructure in relation to the fatigue exposure, without detailing how these changes come about.

Some metal life models approach the fatigue problem from a purely mechanical perspective and treats the metal in a highly idealised way. These models therefore exhibit limitations as regards the description of influences from the operating environment on the characteristics of the bearing material.

A fundamental aspect of the invention is the understanding that during fatigue, the microstructure of materials experiences continuous changes until failure.

According to the invention, by introducing a better material description, further improved predictability power of models can be achieved. Such improved models can be used to better dimension metal objects, e.g. bearings, for adverse operating conditions like environmental influences or combined high load and temperature leading to marginal lubrication conditions. By including additional influencing operating parameters like temperature, internal stress state and speed, improved predictability of e.g. life of mechanical objects, is provided.

According to the invention, fatigue damage is seen as a cumulative small-scale plastic deformation process, being controlled by a thermally activated dislocation climb process. The damage induced is a result of a secondary creep-like dislocation process, where iron self-diffusion controlled climb constitutes the rate controlling step, while the major part of the damage induced is a result of dislocation glide, once the dislocations are freed from the obstacles via the climb process (climb+glide). The damage process is therefore driven by the applied shear stress field and is rate controlled via diffusion-controlled climb.

The present invention is based on a metal physics description of the behaviour of hardened (martensite or bainite) steels. Hardened steels can, from a metal physics point of view, be characterised as being non-equilibrium steel (at equilibrium all steels are soft). Such materials, therefore, behave differently than softer steels. Since fatigue is a result of accumulated damage, induced by incremental micro-plastic deformation in each load cycle, the key to improved predictability lies in understanding the micro-plastic behaviour of hard steels.

While plastic deformation in softer steels is controlled by one deformation mechanism, dislocation glide, or micro-yielding, plastic deformation in hard steels, in their non-equilibrium state, may be induced by two different mechanisms, dislocation glide and dislocation climb. The first deformation mechanism is only active above a given threshold stress level (the micro-yield limit), while the latter one is, although strongly stress dependent, active at all stress levels, not excluding the existence of a lower threshold stress level below which the climb mechanism does not occur. The dislocation climb mechanism is governed by diffusion processes and the fatigue response, therefore, also becomes strongly influenced by temperature, time and internal (hoop and residual) stresses and will therefore also be influenced by the running speed.

There are two classes of residual stresses. The ones that are induced during the manufacturing process and that thus are present from the start of the fatigue process (initial residual stresses), and residual stresses that develop (evolve) in the surface and/or subsurface material during the RCF process as a result of different damage processes. The cumulative plastic strain damage is likely to dominate. Examples of contributing mechanisms to the residual stress evolution are:

-   -   Shakedown, induced by either classic yielding (dislocation         glide) or by a primary creep like deformation process         (dislocation climb+glide). This damage occurs during the early         part of the fatigue process (typically less than one percent of         the fatigue life).     -   Transformation of retained austenite (when applicable) to         martensite or bainite like products (volume increase).     -   Volume changes of the matrix due to accumulation of crystal         defects and due to microstructural changes, followed by:     -   Continuous cumulative plastic strain damage via the         here-discussed (secondary) creep like fatigue damage process         (climb+glide).     -   Microstructural recovery, leading to a partial relaxation of the         residual stress field.

The residual stresses must (just like the load-induced and damage generating dynamic (cyclic) contact stress field) be described as a three-dimensional stress tensor that also varies over time, that is, it is also dynamic (evolving), but on a longer time-scale than the contact-induced stress field. It should also be clear that the residual stress field (tensor), in combination with the hoop stress field (tensor) resulting from interference fitting of inner and outer rings of bearings, or from centripetal forces induced by rotating components, and other superimposed structural stress components. All these various stress fields of different origins, including also the dynamic stress field from the actions of the rolling contact, contribute to the overall, highly dynamic, stress field tensor. From this stress tensor a deviatoric stress field and a normal stress field can be derived. Here the deviatoric stress field drives the fatigue damage process, while the normal stress field influences the rate by which the diffusion-controlled fatigue damage process precedes.

A foundation of the present invention is that the problem of fatigue and factors underlying fatigue are not classic plastic yielding (that is governed by dislocation glide), but a low temperature creep mechanism (that is governed by dislocation climb, itself being rate controlled by iron self-diffusion), rules the fatigue damage process. This latter mechanism is inherently (like all other kinetically controlled chemical and metallurgical processes) temperature dependent and time dependent.

Fatigue exposed components made of hardened and low-temperature tempered steels typically fail due to crack initiation from various defects. In respect to rolling contact fatigue, these defects are present either at the raceway surface (indentations from handling damage or mounting, denting from contaminated running or surface damage induced from improper lubrication, etc.) or within the material, that is, subsurface defects like non-metallic inclusions, pores, pre-existing crack, etc. Defects like these act as local stress raisers leading to a locally accelerated fatigue damage. The shear stress driven and thermally activated fatigue damage processes discussed in this document therefore occur at a higher rate in the steel matrix adjacent to the defects. This results in that the material, at some point in time, will fail locally at these highly stressed points, leading to crack initiation. From this point in time, the fatigue damage process will change to become one of crack propagation till failure (spalling).

In FIG. 1 a device 1 suitable for executing the methods according to present invention is given. It comprises a processor 3 for executing the methods, and input/output means 5, such as a mouse or a keyboard. In an embodiment, it comprises data communication capabilities 7 for receiving and transmitting metal object data. It may also comprise a screen 9 and/or a printer 11 for outputting results from the execution of the methods according to the present invention.

The fatigue test described here is run at high contact pressures, higher than those usually used in bearing applications. The rationale behind this is that the resulting subsurface plastic deformation, which can be observed as an evolving groove at the raceway surface, represents the damage that occur locally at stress raising defects like surface dents or subsurface defects, such as non-metallic inclusions, under typical application conditions.

The surface groove formation is the result of complex processes, where the micro-plastic behaviour is likely to play a central role, but where also other deformation mechanisms and the geometry changes of the contact, induced by the subsurface deformation, etc., also contribute. The grooving model described here describes two contributing parts: The shakedown effect that occurs early in the rolling contact fatigue exposure and which will here be assumed to be cycle independent (immediate); and: The cycle dependent, gradually evolving, grooving mechanism resulting from the subsurface micro-plasticity processes that are described here by the creep-like fatigue damage mechanism.

In terms of shakedown strain damage, the plastic strain, γ_(SD), during the shakedown step is described using power-law relation (the Ludwik relation) γ_(SD) =K< T−T _(u)>^(m) where K is a material constant, T is the shear stress amplitude, τ_(u) is a threshold shear stress level that must be exceeded to induce plasticity and m is the shear stress exponent for the shakedown mechanism. The arrow bracket indicates that if the argument is negative the result is zero.

In terms of steady-state strain damage, the cyclic plastic damage, γ_(N), is modelled as follows (under constant f, temperature and shear stress amplitude). γ_(N) =A _(N) <T−T _(u)>^(c) e ^((−(Q−ΔVσ)/RT)) N/f, where

A_(N) is a material constant, c is the shear stress exponent, Q₀ is the activation energy for the diffusion controlled micro-plastic strain mechanism in absence of internal stress, ΔV is a material constant, σ is the hydrostatic stress state resulting from acting normal internal stresses, N is the number of load cycles (or revolutions), f is the loading frequency (or rotational speed), R is the universal gas constant and T is the local absolute operating temperature. Here, the quotient N/f equals the total running time (t) that is proportional to the “cumulative time under stress” during the rolling contact fatigue exposure.

The core of the model constitutes a “state variable”, that is, a given value of such a variable corresponds to a given micro-plastic damage, independent on which combination of τ, T, t (or N/f) and σ, used to reach the value of this state variable. The state variable may be written: P _(N) =<p ₀ −p _(u)>^(c) e ^((−(Q−ΔVσ)/RT)) N/f, where the shear stress parameters have been replaced by the maximum Hertzian contact pressure (p₀) and a threshold contact pressure (p_(u)), respectively.

Below, two examples of alternative grooving models are shown. The cumulative strain damage for the shakedown and the cycle dependent damage, γ_(sd) and γ_(N), respectively, are used as arguments in a power-law relation: G=(a′·γ _(sd) +b′·γ _(N))^(e′)

An alternative relation is given below, where the difference in character between the shakedown and steady-state strain damage processes is acknowledged by splitting the relation into a direct dependence of the cycle independent shakedown strain and a power-law relation describing the influence of the steady-state strain damage: G=a″·γ _(sd) +b″·γ _(N) ^(e″)

In these relations, the first term accounts for the cycle independent groove formation and the second term accounts for the time or load cycle dependent groove development (sometimes called “ratchetting”). The constants a′, b′ and e′, and a″, b″ and e″, respectively, are regression constants and where the exponents account for the expected non-linearity between the groove depth and the plastic strain due to the groove geometry changes (changing contact geometry, etc.) and other less well defined influences.

The equations above give two different forms of the grooving model. Both these (and other relations) have been assessed in the development of the procedure for the evaluation of groove data. On fitting the experimental data to the selected grooving relation a somewhat modified form of the second equation is used from practical considerations: G=a·

p ₀ −p _(u)

^(m) +b·P _(N) ^(e) where also the shear stress parameters in the shakedown damage relation have been replaced by the maximum Hertzian contact pressure and a threshold contact pressure, respectively, and where the constants a, b and e are used as regression constants. The physically based material parameters (p_(u), m, Q₀, ΔV and c), together with the regression constants, can be evaluated from the experimental information using multi-variable least-square regression analysis.

The grooving model allows also separating the initial (shakedown) damage, and the cyclic or fatigue induced damage. These parameters give additional valuable information on the behaviours of different material variants as basis for decisions on selecting material and processing routes for given applications.

A rig 120 usable when determining values of those constants is schematically shown in FIG. 6. The rig 120 allows running idealised rolling contact fatigue (RCF) tests where essentially only the steel matrix response to the RCF conditions are evaluated. In an embodiment, the specimen is a cylindrical roller 122 (diameter 18 mm and length 30 mm), which is exposed to near point contacts by two ceramic (Si₃N₄) balls 124 (diameter 19 mm), running along the same track on the roller circumference, giving two stress cycles per specimen revolution. Also, lubrication means 126 for lubricating the near point contacts is also provided. Three bearings 126 are provided to support the balls 124 and a non-driven wheel 128 supports the specimen 122.

This allows controlling closely the operating conditions, allowing varying contact pressure, varying temperature and varying rotational speed. Full film lubrication is maintained at all times. The fatigue response of the steel matrix can be evaluated using different techniques like metallography (optical and electron microscopy) and X-ray techniques, etc., methods that can be characterised as “indirect”. In the method described here a more direct method of determining the cumulated plastic damage is applied, namely recording the groove formation at the raceway surface, resulting from the subsurface small-scale plastic damage during the RCF exposure.

A typical test programme for a material of interest consists of a test matrix involving a number of experiments usually run at constant speed and at varying contact pressure levels (at least two) and varying temperature (at least two levels). This allows evaluating the basic material parameters, or constants, including p_(u), m, Q₀, and c. If also the material constant ΔV should be evaluated, additional experiments with varying internal stress levels such as (tensile) hoop stress (for example obtained by use of a compound roller with a shrink fitted shell on a roller core) or (usually compressive) residual stress (introduced during manufacturing of the specimen), have to be run (at least two internal stress levels have to be introduced). For each of the selected test conditions the evolving surface groove is evaluated using surface profiling techniques and a defined profile parameter (for example, groove shoulder to groove bottom or groove depth below the original surface) is defined to represent the groove depth. This parameter is measured either in-situ, using on-line interferometry, or intermittently for different exposure times, ranging from short (seconds or minutes) to long (hours or days).

In FIG. 3, a diagram indicating the relation between subsurface plastic fatigue damage observed in the form of a surface groove, plotted versus number of load cycles, and temperature, is given. This is an example of the correlation between the grooving model and experimental data from the testing of martensitically induction hardened, unalloyed, medium carbon bearing steel. The following set of material constants were found, as a result of the simultaneous multi-variable regression analysis of all data points (n=24):

m=2.5

p_(u)=580 MPa

c=16.3

Q₀=100 500 J/mol

In this case the ΔV value could not be evaluated since all tests were run without changing the internal stress parameter. The groove model fitting parameters took here the values: a=2.44·10⁻¹⁰, b=1.48·10⁻¹³ and e=0.254, with an overall correlation coefficient r²=0.986. The diagram shows the grooving response after testing at two levels of contact pressure, 3.2 and 4.9 GPa, and at two temperatures, 65 and 100° C. Markers represent measured groove depths and the curves are the result of a simultaneous, multi-variable regression analysis of all data points shown (n=24, r²=0.986). Also applicable to other figures, (Exp) denotes experimental data and (Mod) denotes data as originated from the grooving model.

In FIGS. 4A and 4B, diagrams indicating the relation between subsurface plastic fatigue damage observed in the form of a surface groove, plotted versus number of load cycles, and temperature, are given. This is an example of the correlation between the grooving model and experimental data from the testing of martensitically through hardened, low-alloyed, high-carbon bearing steel. The following set of material constants were found as a result of the simultaneous multi-variable regression analysis of all data points (n=40):

-   -   m=10     -   p_(u)=1200 MPa     -   c=13.6     -   Q₀=135 000 J/mol

In this case the ΔV value could not be evaluated since all tests were run without changing the internal stress parameter. The groove model fitting parameters took here the values: a=4.13·10⁻³⁷, b=9.55·10⁻¹⁰ and e=0.258, with an overall correlation coefficient r²=0.979. The diagrams show the grooving response after testing at a range of contact pressures and at two temperatures, 65 and 100° C. Markers represent measured groove depths and the curves are the result of a simultaneous, multi-variable regression analysis of all data points shown (n=40, r²=0.979).

In FIG. 5, an example of the correlation between the grooving model and experimental grooving data of martensitically hardened, chromium and vanadium alloyed tool steel of bearing quality is shown. The following set of material constants were found for this secondary hardened tool steel, as a result of the simultaneous multi-variable regression analysis of all data points (n=35):

-   -   m=6     -   p_(u)=2500 MPa     -   c=17.2     -   Q₀=106 000 J/mol

In this case the ΔV value could not be evaluated since all tests were run without changing the internal stress parameter. The groove model fitting parameters took here the values: a=9.34·10⁻²², b=4.16·10⁻⁷ and e=0.132, with a correlation coefficient of r²=0.988. The diagram shows the grooving response after testing at two levels of contact pressure, 4.0 and 5.0 GPa, and at two temperatures, 75 and 100° C. Markers represent measured groove depths and the curves are the result of a simultaneous, multi-variable regression analysis of all data points shown (n=35, r²=0.988). 

The invention claimed is:
 1. A method of indicating a fatigue damage rate of a hardened metal object in relation to load cycles, N, exerted on the hardened metal object, comprising the step of: calculating, using a processor, the fatigue damage rate based on an effective activation energy parameter for a dislocation climb process, Q, a shear stress amplitude, τ, an absolute local temperature of the hardened metal object, T, and a load frequency, f.
 2. The method according to claim 1, wherein the step of calculating includes calculating the fatigue damage rate according to dγ/dN=A<τ−τ _(u)>^(c) e ^((−kQ/RT)) /f, where: A is A₀(H_(Vref)/H_(V))^(d), where A₀ is a constant greater than zero and lesser than 10⁻⁸, H_(Vref) is the Vickers hardness for a reference metal, H_(V) is the Vickers hardness for the hardened metal object, <τ−τ_(u)> is zero for τ≦τ_(u) and τ−τ_(u) for τ>τ_(u), where τ_(u) is the fatigue limit for the hardened metal object, c is a constant in the interval 6 to 22, d is a constant in the interval 6 to 22, k is a constant between 0.50 and 1.50, T is the absolute temperature of the hardened metal object at the point of contact, and R is the universal gas constant.
 3. The method according to claim 1, wherein the step of calculating includes calculating the fatigue damage rate according to dγ/dt=A<τ−τ _(u)>^(c) e ^((−kQ/RT)), where: A is A₀ (H_(Vref)/H_(V))^(d), where A₀ is a constant greater than zero and lesser than 10⁻⁸, H_(Vref) is the Vickers hardness for a reference metal, H_(V) is the Vickers hardness for the hardened metal object, <τ−τ_(u)> is zero for τ≦τ_(u) and τ−τ_(u) for τ>τ_(u), where τ_(u) is the fatigue limit for the hardened metal object, c is a constant in the interval 6 to 22, d is a constant in the interval 6 to 22, k is a constant between 0.50 and 1.50, T is the absolute temperature of the hardened metal object at the point of contact, and R is the universal gas constant.
 4. The method according to claim 1, wherein the temperature, T, is in the interval of zero degrees Celsius to 250 degrees Celsius.
 5. The method according to claim 1, wherein the temperature, T, is in the interval of 15 degrees Celsius to 200 degrees Celsius.
 6. The method according to claim 1, wherein the temperature, T, is lesser than 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.
 7. The method according to claim 1, further comprising the step of calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q ₀ −ΔVσ, where ΔV is a material constant, and σ is a normal stress field.
 8. The method according to claim 7, wherein σ is hydrostatic stress.
 9. The method according to claim 8, wherein σ is calculated according to σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are normal stress components.
 10. A method of indicating fatigue damage, γ, of a hardened metal object at a temperature substantially corresponding to operating conditions of the hardened metal object and under a constant stress amplitude, comprising the step of: calculating, using a processor, the fatigue damage, γ, based on an effective activation energy parameter for a dislocation climb process, Q, a shear stress amplitude, τ, an absolute local temperature of the hardened metal object, T, a load frequency, f, and a number of load cycles, N, exerted on the hardened metal object.
 11. The method according to claim 10, wherein the step of calculating includes calculating γ according to γ=A<τ−τ _(u)>^(c) e ^((−kQ/RT)) N/f, where A is A₀(H_(Vref)/H_(V))^(d), where A₀ is a constant greater than zero and lesser than 10⁻⁸, H_(Vref) is the Vickers hardness for a reference metal, H_(V) is the Vickers hardness for the hardened metal object, <τ−τ_(u)> is zero for τ≦τ_(u) and τ−τ_(u) for τ>τ_(u), where τ_(u) is the fatigue limit for the hardened metal object, c is a constant in the interval 6 to 22, d is a constant in the interval 6 to 22, k is a constant between 0.50 and 1.50, T is the absolute temperature of the hardened metal object at the point of contact, and R is the universal gas constant.
 12. The method according to claim 10, wherein the step of calculating includes calculating γ according to γ=A<τ−τ _(u)>^(c) e ^((−kQ/RT)) t, where: A is A₀(H_(Vref)/H_(V))^(d), where A₀ is a constant greater than zero and lesser than 10⁻⁸, H_(Vref) is the Vickers hardness for a reference metal, H_(V) is the Vickers hardness for the hardened metal object, <τ−τ_(u)> is zero for τ≦τ_(u) and τ−τ_(u) for τ>τ_(u), where τ_(u) is the fatigue limit for the hardened metal object, c is a constant in the interval 6 to 22, d is a constant in the interval 6 to 22, k is a constant between 0.50 and 1.50, T is the absolute temperature of the hardened metal object at the point of contact, and R is the universal gas constant.
 13. The method according to claim 10, wherein the temperature, T, is in the interval of zero degrees Celsius to 250 degrees Celsius.
 14. The method according to claim 10, wherein the temperature, T, is in the interval of 15 degrees Celsius to 200 degrees Celsius.
 15. The method according to claim 10, wherein the temperature, T, is lesser than 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.
 16. The method according to claim 10, further comprising the step of calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q ₀ −ΔVσ, where ΔV is a material constant, and σ is a normal stress field.
 17. The method according to claim 16, wherein σ is hydrostatic stress.
 18. The method according to claim 17, wherein σ is calculated according to σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are normal stress components.
 19. A method of indicating at least one of an effective activation energy parameter for the a dislocation climb process, Q, a shear stress amplitude, τ, an absolute local temperature of a hardened metal object, T, a load frequency, f, and a number of load cycles, N, exerted on the hardened metal object, at a temperature substantially corresponding to operating conditions of the hardened metal object and under a constant stress amplitude, comprising the step of calculating, using a processor: N=Cf/(<τ−τ_(u)>^(c) e ^((−kQ/RT))), where: C=C₀(H_(V)/H_(Vref))^(d), where C₀ is a calibration constant representing operating conditions for a reference metal object having a known fatigue life in terms of estimated maximum number of load cycles, H_(Vref) is the Vickers hardness for a reference metal and H_(V) is the Vickers hardness for the hardened metal object, <τ−τ_(u)> is zero for τ≦τ_(u) and τ−τ_(u) for τ>τ_(u), where τ_(u) is the fatigue limit for the hardened metal object, c is a constant in the interval 6 to 22, d is a constant in the interval 6 to 22, k is a constant between 0.50 and 1.50, T is the absolute temperature of the hardened metal object at the point of contact, and R is the universal gas constant.
 20. The method according to claim 19, wherein the temperature, T, is in the interval of zero degrees Celsius to 250 degrees Celsius.
 21. The method according to claim 19, wherein the temperature, T, is in the interval of 15 degrees Celsius to 200 degrees Celsius.
 22. The method according to claim 19, wherein the temperature, T, is lesser than 0.4T_(m), where T_(m) is the absolute melting temperature of the hardened metal object.
 23. The method according to claim 19, further comprising the step of calculating the effective activation energy parameter for the dislocation climb process, Q, according to Q=Q ₀ −ΔVσ, where ΔV is a material constant, and σ is a normal stress field.
 24. The method according to claim 23, wherein σ is hydrostatic stress.
 25. The method according to claim 24, wherein σ is calculated according to σ=(σ_(xx)+σ_(yy)+σ_(zz))/3, where σ_(xx), σ_(yy), σ_(zz) are normal stress components.
 26. A non-transitory computer-readable medium storing instructions that, when executed by one or more processors of a computing system, are configured to cause the computing system to perform operations comprising: calculating a fatigue damage rate based on an effective activation energy parameter for a dislocation climb process, Q, a shear stress amplitude, τ, an absolute local temperature of a hardened metal object, T, and a load frequency, f so as to indicate fatigue damage rate of the hardened metal object in relation to load cycles, N, exerted on the hardened metal object; or calculating fatigue damage, γ, based on the effective activation energy parameter for the dislocation climb process, Q, the shear stress amplitude, τ, the absolute local temperature of the hardened metal object, T, the load frequency, f, and the number of load cycles, N, exerted on the hardened metal object so as to indicate the fatigue damage,γ, of the hardened metal object at a temperature substantially corresponding to operating conditions of the hardened metal object and under a constant stress amplitude; or both. 